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O ' Clusters of galaxies and the large scale filaments that connect neighboring clusters are ex- 

, pected to be sites of acceleration of charged particles and sources of non-thermal radiation 

■ from radio frequencies to gamma rays. Gamma rays are particularly interesting targets of 

investigation, since they may provide precious information on the nature and efficiency 
of the processes of acceleration and magnetic confinement of hadrons within clusters 
of galaxies. Here we review the status of viable scenarios that lead to the production 
i of gamma rays from large scale structures and are compatible with the multifrequency 

J.^ ' observations that are already available. We also discuss the possibility of detection of 

gamma rays with space-borne telescopes such as GLAST and ground based Cherenkov 
telescopes, and the physical information that may be gathered from such observations. 



1. Introduction 

Clusters of galaxies are large collections of galaxies, gas and especially dark mat- 
ter, confined in a volume of a few Mpc 3 . Clusters and the filaments that connect 
them represent the largest structures in the present universe in which the gravita- 
tional force due to the matter overdensity overcomes the expansion of the universe. 
This condition allows clusters to collapse and eventually go through a process of 
virialization at the present cosmic time. 

Rich clusters have typical total masses of the order of 1O 15 M , mostly in the 
form of dark matter, while ~ 1% is in the form of galaxies and ~ 5% is in the form 
of a hot (T - 10 8 K), tenuous (n gas - 10~ 3 - 10~ 4 cm~ 3 ), X-ray emitting gas^. 
In terms of energy density, the gas in the cluster is typically heated to roughly the 
virial temperature, but there is room to accomodate, in principle, a non negligible 
amount of non-thermal energy, in the form of accelerated particles. 
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A non-thermal component is in fact observed from several clusters of galaxies, 
mainly in the form of a tenuous radio emission^. In some cases a hard (non-thermal) 
X-ray emission is also observed^, whi le more debate exists on the non-thermal origin 
of a soft X-ray and UV emission The radio emission has a steep spectrum and 
appears mainly in two forms, halos and relics, the former being very extended, 
fairly regular and at the cluster center, and the latter being typically elongated and 
located in the outskirts of clusters. 

Roughly 20% of clusters with luminosity in the 0.1-2.4 keV band > 5 x 10 44 
erg/s are observed to have radio halos In general radio halos are rather rare 
phenomena^, which means that their typical lifetime is relatively short (~ 1 Gyr) 
- . Such radio emission can only be interpreted as the result of synchrotron emission 
in the intracluster medium, thereby confirming that 1) relativistic electrons are 
present, 2) the intracluster medium (ICM) is magnetized and 3) the magnetic field 
structure is topologically complex. The third point is inferred from the fact that the 
synchrotron emitting particles should be diffusively trapped on the scale of radio 
halos 

Independent evidence for turbulent magnetic fields in the ICM comes from Fara- 
day rotation measurements (RM) from radio sources located inside or behind clus- 
|- erg | ll | 12 l This method leads to values of the field of the order of several /xG^^, but 
the results are affected by uncertainties in the topology of the field and the spatial 
distribution of the electron gas, as well as by the subtraction of the intrinsic RM at 
the source. Smaller fields, of the order of a few tenths of /j,G, are obtained from the 
combination of radio and hard X-ray measurements, in the few cases in which data 
at these wavelengths are available from the same cluster. This method relies on the 
assumption that the diffuse radio emission and the hard X-ray excess are cospa- 
tial and produced by the same population of relativistic electrons via synchrotron 
emission and inverse Compton scattering respectively Although it is generally 
accepted that at least a fraction of the ICM is magnetized at /iG level, the intrinsic 
uncertainties in both these methods do not allow yet a precise determination of the 
magnetic field and its spatial profile 

The electrons responsible for the radio emission may be accelerated in a variety of 
ways and in a variety of sources within the ICM. However, it is customary to classify 
the models of the origin of the non-thermal activity in primary electron models, 
secondary electron models and reacceleration models, depending on whether the 
radiating electrons are 1) accelerated in specific sites (e.g. shocks, AGNs, galaxies), 
2) produced as secondary products of hadronic interactions, or 3) continuously 
reaccelerated from a pre-existing population of non-thermal seeds in the ICM. 

Shock waves associated to the process of structure formation are expected to be 
the sites where most cosmic rays in clusters arc accelerated. In particular, the large, 
strong shocks that bound the motion of matter in filaments and during accretion 
onto clusters are believed to be very efficient accelerators. 

To date, no non-thermal activity has reliably been observed from accretion 
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shocks or large scale filaments @. As we discuss later, these structures may how- 
ever be detected in the near future or lead to phenomena that may be observable 
from within the cluster volume. 

No cluster of galaxies has been firmly detected in gamma rays so far S3. De- 
spite this non detection, two ideas contribute to reinforce the beli ef t hat clus ters of 
galaxies can be interesting as high energy radiation emitters: in and it was 
first understood that the bulk of cosmic rays accelerated within the cluster volume 
would be confined there for cosmological times, thereby enhancing the possibility of 
inelastic proton-proton collisions and consequent gamma ray production through 
the decay of neutral pions. 

In I2£U it was recognized that large scale shocks associated with the formation of 
structures in the universe, may accelerate electrons to TeV energies, implying that 
high energy emission would occur due to the upscattering of the photons of the 
cosmic microwave background to gamma ray energies, through inverse Compton 
scattering (ICS). The detection of gamma rays from clusters of galaxies, or more in 
general from large scale structures would represent an important and crucial step 
forward towards a better understanding of the processes of acceleration of particles 
in the ICM and cosmic ray confinement. 

The measurement of the gamma ray flux associated with clusters can allow us 
to discriminate among vastly different scenarios for the origin of the non-thermal 
activity associated with the formation of large scale structures in the universe. 

This review is organized as follows: in Sj5]we discuss the important phenomenon 
of cosmic ray confinement in the ICM and the main energy loss processes for leptons 
and hadrons. In fJ3]wc summarize our understanding of the sources and acceleration 
sites for cosmic rays in clusters of galaxies. 5E1 is devoted to a discussion of the 
formation of large scale shocks during mergers or accretion and to assessing their 
role as cosmic ray accelerators. In £[5]we discuss some recent developments in the 
theory of particle acceleration at shock waves, and their possible role for large scale 
structure formation shocks. In Sj6]we summarize the main predictions on the fluxes 
of gamma rays from clusters of galaxies, both in the GeV and TeV energy domains. 
The contribution of clusters of galaxies to the diffuse gamma ray background is 
discussed in [J7l Our conclusions are provided in EjSj 

2. Dynamics of accelerated particles in clusters of galaxies 

The propagation of charged particles (nuclei and electrons) injected in the ICM is 
mainly determined by diffusion, convection and energy losses. 

As long as the turbulent part of the magnetic field remains smaller than the 
large scale field, the random walk in pitch angle that the particles perform, and the 

a Recently, a de tection of radio waves from a cosmic shock surrounding the cluster A3376 has been 
claimed in ED. However, alternative explanations, for instance associated with outgoing merger 
shock fronts may be equivalently consistent with the data. 
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spatial diffusion that follows are described by the so-called quasi-linear theory^^. 
The diffusion coefficient for scattering of particles with AlfVen waves is therefore 

r L (p)v(p) 1 B 2 

where p is the particle momentum, r^(p) = pc/ZqB is the Larmor radius for a par- 
ticle of charge Z in a magnetic field B, v(p) ~ c is the particle velocity and J 7 (k(p)) 
is the energy density in the form of turbulent field on a scale k that interacts reso- 
nantly with particles with momentum p oc 1/k. We introduced the power spectrum 
of the perturbations of the magnetic field P(k), such that J fc °° dkP(k) = £_B 2 and 
£ < 1. The minimum wavenumber k m i n — 2Tr/L max is related to the largest scale 
in the turbulent magnetic field. In clusters it is usually believed that L max ~ 100 
kpc. For a Kolmogorov spectrum P(k) oc fc~ 5 / 3 , the diffusion coefficient reads 

D(E) « 7 x l^EiGeVylz^B-^L 2 ^ cmV\ (2) 

where is the magnetic field in units of /iG and iioo = ( lookpc ) ■ F& ra day rotation 
measures (RM) support a scenario in which most of the magnetic fi eld energy is 
associated with a turbulent field tangled on both small and large scales At small 
scales (few kpc) a Bayesian analysis of the RM suggests that the power spectrum 
of the field follows a Kol mog orov scaling 1221, although the situation at larger scales 
might be more complex ■■ ~. In the absence of better estimates or measurements 
of the intracluster magnetic field strength and topology, here we assume £ ~ 1. 
Larger values of £ are hardly treatable within quasi-linear theory but qualitatively 
they lead to more effective confinement, although it is not even clear whether the 
particle propagation remains purely diffusive in such a limit. 

The diffusion time on the scale i? ~ 1 Mpc of the whole cluster is t^// = 
R 2 /4D(E), which exceeds the age of the cluster (to ~ 10 Gyr) for energies E < 
E c = 1 BfiL^Q TeV. Particles with energies lower than E c are confined within the 
cluster volume, a phenomenon which is unique for clusters of galaxies and is mainly 
due to their very large physical size I 18 | 19 | 

The prediction of diffusive confinement of cosmic rays within the cluster volume 
is very strongly dependent upon the choice of the diffusion coefficient and even the 
specific assumptions on the parameters involved, due to both the energy content 
in the form of turbulent magnetic field (the parameter £) and the spectrum of 
fluctuations in the B-field. For instance a spectrum with P(k) oc fc -3 / 2 would lead 
to E c = 4.7 x 10 5 £ B^L^qq GeV. The dependence on the coherence scale of the field 
is also important, the confinement being more efficient for smaller values of L max . 

While diffusing in the intracluster magnetic fields, charged particles are subject 
to energy losses: for electrons the dominant loss channels are ICS off the photons 
of the cosmic microwave background (CMB), synchrotron emission, and Coulomb 
losses (at low energies) . Protons on the other hand lose energy mainly through pp 
inelastic scattering when their kinetic energy is higher than ~ 300 MeV. At lower 



February 5, 2008 15:34 



Gamma rays from clusters of galaxies 5 

energies Coulomb losses become important. The time scales for losses due to the 
combination of these processes are illustrated in Fig. [TJ The curve in the middle 
is the time scale of losses for protons, while the bottom curves represent the loss 
time scale for electrons in a mean magnetic field B = 1/iG (solid line) and 3/iG 
(dashed line). In the same plot we also show the confinement time for the case of 
a Kolmogorov spectrum with £ = 1 and L max = 100 kpc and mean magnetic field 
B = 1/iG (solid line) and 3/iG (dashed line). 

One consequence of the confinement of the hadronic component in the ICM 
is that the grammage traversed by cosmic rays is increased to the level that a 
substantial gamma ray production can in principle take place. The grammage for 
a cosmic ray proton with energy £ in a cluster is X(E) — X c — m p n gas cto wl- 
10 g cm -2 for E < E c (here m p is the proton mass and n gas is the gas density) . For 
E > E c the grammage decreases with energy. The value of X c should be compared 
with the nuclear grammage X nuc w 50 g cm -2 , which implies that an appreciable 
amount of the energy of a low energy cosmic ra y is los t in inelastic pp scatterings, 
leading to gamma ray and neutrino production I 18 | 19 | anc j ^ ^ ne generation of a 
secondary population of electrons and positrons. 




Log(pc/GeV) 

Fig. 1. Time scale for energy losses of protons (middle) and electrons (bottom). The lines in 
the top part are the diffusion time scales for protons for a Kolmogorov spectrum of magnetic 
fluctuations. The magnetic field is 1/iG (solid) and 3fiG (dashed line). The feature at 1—3 GeV in 
the time scale for energy losses of protons is due to the combination of the shapes of the inclusive 
pp cross section for 7r + and ir~ . 



From Fig. [T] we can also conclude that electrons accelerated at specific sites can 
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hardly radiate their energy far from the production site 1=*h This results in the fact 
that the region of radio emission is very localized, quite unlike what is observed 
in radio halos. Radio emission at ~ 1 GHz is typically generated by electrons with 
energy ss 10B^ ^ 2 GeV, and at this energy the loss time is « 10 8 years. The distance 
covered by diffusion during the loss time is tiny compared with the observed size 
of the radio emitting region (several hundred kpc at this frequency) . The electrons 
could move farther from the production site due to advection (for instance behind 
a shock), but this would need that fast shocks, with Mach number larger than 
4-5, should cross the cluster cores, and still the radio emission would decay after 
< 100 million years, which would make radio halos less frequent than observed ^Sl. 
These arguments lead to the requirement that electrons are generated or accelerated 
everywhere in the cluster, which leaves us with two possibilities: eithe r the r adiating 
electrons have a secondary origin, being produced by pp scatterings l ^7 | 28 ] or they 
are continu ously reaccelerated in situ through a second order Fermi mechanism 
| 29 | 30 | 31 110| ;g e i ow we discuss these possibilities in some more detail. 



3. Sources of cosmic rays and models of particle acceleration 

Whether the radiating electrons are accelerated at specific sites, reaccelerated from 

pre-existing seeds or produced as secondary products of hadronic interactions, a 

class of sources of accelerated particles (electrons, hadrons or both) must exist in the 

ICM. Particle acceleration can take place in several places in the ICM, from ordinary 

galaxies, which contribute a total luminosity in the form of protons ~ 3x 10 42 erg s^ 1 , 

to active g alaxies, which m ay plausibly contribute ~ 10 45 erg s" 1 over periods of 
„ 10 8 ycar J19l32l28l33l34l 

The process of structure formation contributes the most to the ener- 
getics of clusters, both i n the thermal EH an d non-thermal components 
| 18 | 19 | 36 | 37 | 42|4 3 38 39 40|41] Mergers between two or more clu sters are observed 
to heat clusters I 44 | 45 | anc j p 0S gi D iy accelerate particles I 46 | 47 | 48 l^ w hii e ^ j s harder 
to measure the effect of the continuous infall of material onto clusters from the 
colder outskirts. The two processes are not really independent, but for simplicity 
we will refer to the former as mergers and to the latter as accretion. Mergers provide 
the largest contribution to the heating of the ICM, mainly through the production 
of weak shocks that transform an appreciable fraction of the relative kinetic energy 
of the merging clusters into thermal energy. Accretion is less important in this re- 
spect but it is likely to play an important role for non-thermal phenomena (see next 
section). 

The total energy in a merger of two clusters with roughly equal mass M = 
10 14 M Q is E w 10 63 erg. A few of these events may lead the resulting cluster to reach 
the virial temperature or close to it. It is important to realize that the majority of 
the mass relevant for this estimate is in the form of dark matter, while only baryons 
(and electrons) are heated. In the assumption that the gaseous components of the 
initial clusters are at the virial temperature, it is easy to show that the merger occurs 
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at supersonic relative speed, therefore implying the formation of shock waves. The 
same conclusion can be reached for the accreting material^. If a fraction of this 
energy can be converted at the shocks into non-thermal particles through a first 
order Fermi process, then the ICM may be populated with a large amount of non- 
thermal particles, which may potentially become important even for the dynamics 
of the cluster ^1 Both electrons and protons can be accelerated at shocks, and 
secondary electrons can in turn be produced during inelastic collisions of cosmic 
ray protons in the ICM. 

In terms of the morphology of the resulting non-thermal emission these scenarios 
of injection of relativistic electrons are not equivalent. If electrons are accelerated at 
merger shocks, the resulting radio emission is found to have a filamentary structure 
tracing the position of the shock ^01, rather than the regular, symmetric structure 
observed in radio halos. That is because, as discussed in Sec. [21 electrons cannot 
diffuse far away from their acceleration site, due to the short cooling ti me. Radio 
relics are indeed interpreted as radio emission triggered by merger shocks I 46 | 51 | 52 ] 

Extended and fairly regular radio emission is predicted when the radiating 
electrons ar e of secondary origin, since parent protons can diffuse on large scales 
| 27 | 28 | 53 [50| However two aspects of radio halos are difficult to explain in this sce- 
nario: 1) since all clusters have suffered mergers (hierarchical scenario) and hadrons 
are mostly confined within clusters, extended radio emission should be basically 
observed in almost all clusters. Radio n .^P. s . _ ar e presently detected only in a fraction 
of massive and merging clusters |7 | 8 | 54 | 55 | 2 | 3 [ 2) Some of the observed r adio halos 
are found to have a synchrotron spectrum with a cutoff at a few GHz | 56 | 57 | 58 | 
corresponding to a cutoff in the electron spectrum. This, together with the fact 
that radio halos are very extended, would require art ificial a nd unrealistic assump- 
tions in order to be consistent with secondary models l ^ l ^^ l. On the other hand the 
observational facts illustrated above are mainly based on a handful of clusters and 
might reveal themselves as incomplete when a more significant sample will become 
available. It is possible that some halos, possibly not yet identified, may have sec- 
ondary origin, despite the fact that the ones that we are now aware of can hardly 
be exlained in this way. 

Gamma ray emission from clusters of galaxies is still believed to be the most 
effective way to detect the presence and confinement of cosmic ray hadrons in the 
intracluster medium. 

Present observations of extended radio emission and, when present, of hard X- 
ray emission, are best explained by models of continuous reacceleration of electrons 
due to particle- wave interactions. A recent review can be found in 1221. Here we limit 
ourselves to provide a summary of the basic aspects of this class of models. 

Re acceleration models are basically models of second order Fermi acceleration 
| 60 | 61 |^ m w j 1 j c j 1 charged particles are accelerated stochastically due to the random in- 
teraction of the particles with perturbations (waves) in the structure of the magnetic 
field. These wave-particle interactions typically require that mildly non-thermal par- 
ticles (lorentz factors 7 ~ 100 — 500) are already present in the region (seed parti- 
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cles). In the case of clusters of galaxies these seeds could be provided by the past 
activity of ac tive galaxies in the ICM or more plausibly by the past merger history 
of the cluster EH. Alternatively the seeds could be secondary products of hadronic 
interactions ^2l. 

As stressed above, the extended radio emission from clusters of galaxies is a 
rather rare phenomenon, hinting to the fact that the episodes of radio activity are 
relatively short (~ 1 Gyr). In the case of wave-particle interactions, this scale is 
automatically recovered as a combination of the time necessary for the cascading of 
turbulence from large scales (300-500 kpc) to the small scales relevant for particle 
acceleration, and of the cluster-cluster crossing time. 

Most scenarios of particle reacceleration in the ICM are based on particle in- 
teractions with either large scale compressible modes (magnetosonic) or small scale 
Alfven modes. Resonant particle acceleration due to Alfven modes is one of the 
most common mechanisms used in particle acceleration in astrophysics, as the res- 
onant coupling with fast particles is the main channel of energy dissipation of these 
modes. Alfvenic reacceleration has been used to explain the origin o f the e mitting 
electrons in radio halos l ®3 | 64 [ Detailed time dependent calculations I 65 | 62 | SU gg es t 
that efficient Alfvenic reacceleration of relativistic electrons might operate in the 
ICM provided that the energy budget of relativistic hadrons is below 3-5% of that of 
the thermal ICM and that a few percent of this thermal energy is channelled into the 
modes during a lifetime of a radio halo. Relativistic hadrons affect the reacceleration 
of fast electrons because they provide the dominant contribution to the damping of 
Alfven modes in the ICM: the acceleration of such hadrons increases the damping 
of the mo des and this limits the acceleration of relativistic electrons (wave-proton 
boiler) ^^51. The injection process of Alfven modes in the ICM and the assumption 
of isotropy of these modes represent the most important source of uncertainties in 
this class of models. Alfven waves couple with relativistic electrons (and protons) 
on very small spatial scales I ~ 2np/ (Qm) (p, m and Q being the momentum, mass 
and Larmor frequency of particles) and this requires a surprisingly efficient Alfvenic 
cascade if these modes arc connected with the merger-driven turbulence on larger 
scales. In addition, the cascading of MHD turbulence in nature usually develops 
into a n anisotropic spectrum on scales significantly smaller than the injection scale 
| 66 | 67 | anc j ^his should strongly reduce the efficiency of particle acceleration EH. On 
the other hand, injection of Alfven modes at (quasi-) resonant scales might occur as 
a consequence of the cascading of fluid turbulence or of several kinds of plasma 
instabilities the physical details and efficiency of these mechanisms are however 
yet to be understood. 

Some of the problems encountered with Alfven modes may be solved if the 
interaction of particles with compressible modes is considered, although new is- 
sues are also raised. Compressible modes do not require injection on small spatial 
scales: a substantial fraction of the cluster turbulence is expected to be in the 
form of large scale compressible isotropic modes, which might efficiently acceler- 
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ate fast particles via r esonant Transit Time Damping (TTD) and non-resonant 
turbulent-compression | 41 | 71 | 72 ] rp^g ma j n source f uncertainty in these models 
is our ignorance of the viscosity in the ICM. Such viscosity could in fact severely 
damp compressible modes on large scales, inhibiting particle acceleration. On the 
other hand, turbulence in the ICM is subsonic but strongly super- Alfvenic, and 
in these conditions the bending of the magnetic field lines is expected to main- 
tain the effective i on-i on mean free path much smaller than that in the classical 
unmagnetized case ^1. This phenomenon should limit the viscous dissipation, so 
that collisionless dampings with the thermal and relativistic particles should repre- 
sent the main source of turbulent dissipation. In this case detailed calculations of 
resonant -TTD and non-resonant stochastic coupling with fast modes in the ICM 
suggest that the resulting efficiency of the particles re-acceleration could be suf- 
ficient to power radio halos, provided that the rms velocity of these compressible 
eddies is V^ ms ~ (0.1 — 0.3)Cg ^1 Here, at variance with the Alfvenic case, the 
acceleration efficiency of electrons and positrons is not affected by the presence of 
relativistic hadrons since the damping of compressible modes in the I CM is lar gely 
dominated by the TTD-coupling with thermal electrons and protons [ 74 | 72 | 41 | 

In models with stochastic reacceleration, the spectrum of electrons has typically 
a cutoff at Lorentz factors < 10 5 , therefore no appreciable ICS gamma rays may 
be expected. However, as pointed out in E2 : , the role of protons in these models 
may be very important and one may reasonably expect gamma ray emission from 
either reaccelerated protons or through ICS of the secondary electrons generated in 
inelastic pp collisions. Moreover, stochastic reac celerat ion may significantly modify 
the shape of the spectrum of relativistic hadrons | 65 | 72 | thereby changing the shape 
of the spectrum of neutral pions and gamma rays. Further investigations of these 
effects are certainly needed and desirable. 



4. Large scale shocks during cluster formation 

As discussed above, most of the non-thermal energy content of clusters is believed to 
be associated with the process of large scale structure formation. We briefly review 
this process here, with particular emphasis on the issues which are most important 
for particle acceleration and related non-thermal phenomena in clusters. 

A wide consensus has been reached in the community on the hierarchical model 
of structure formation, which predicts that increasingly larger halos are formed by 
merging of smaller halos. On average smaller structures are therefore formed at 
earlier times than larger structures. During mergers between clusters, shocks are 
driven by gravity in the diffuse baryonic component, which is heated up to X-ray 
temperatures^. High resolution X-ray observations of clusters of galaxies performed 
with Chandra and XMM have provided us with direct eviden ce for disturbed X-ray 
morphology and hot shocks associated with merging clusters I 75 | 76 | 

While a complete understanding of the process of structure formation can only 
be achieved through numerical N-body simulations, a simple analytical description 
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is useful in that it allows us to check the main numerical results versus some basic 
physical expectations. These analytical descriptions come in different flavors and are 
widely discussed in the literature. Historically, t he fi rst approach to the proble m was 
proposed by Press & Schechter (hereafter PS) and further developed in l ?8 | 79 | 
among others. 

The PS formalism provides us with the differential comoving number density 
of clusters with mass M at cosmic time t, n(M,t), and with the rate at which 
clusters of mass M merge at a given time t to give a cluster with final mass M' , 
R{M, M', t). These two distributions can be randomly sampled to construct Monte 
Carlo realizations of the merger history of a cluster. This procedure is illustrated in 
the left panel of Fig. 2 ESI f or a cluster with present mass of 10 15 M Q . The history 
is followed backward in time to a redshift z = 3. The big jumps in the cluster 
mass correspond to major merger events, in which two subclusters collide to form 
a bigger object, while smaller jumps correspond to accretion events in which the 
cluster mass is increased by a small amount. 

The PS formalism has been extensively used in th e past to inve stigate the heating 
of the ICM and other thermal properties of clusters EmU . In EH such approach has 
been generalized in order to describe the acceleration of particles at merger shocks 
and to investigate the implications for clusters non-thermal emission. The novelty 
in the approach resides in the combination of the PS formalism with a recipe to 
evaluate the Mach number of merger related shocks. An accurate determination of 
the shock strength is of crucial importance for predicting the non-thermal behavior 
of clusters of galaxies, since the spectrum of the particles accelerated at a shock 
is, at least in the test particle approximation, uniquely determined by its Mach 
number. 

The Mach number of merger-related s hoc ks c an be estimated by using an ap- 
proach introduced in its original version in EH and ^Sl. The basic assumption is that 
the two merging clusters, with masses M\ and Mi are virialized to start with. The 
virial radius for each cluster can be written as follows: 

(GM \ ^ ( Mi \ 

ioon,W+ */..)■ ) ~ % \wna) t^'f'r'Mpc, (3) 

where i = 1,2, Hq is the Hubble constant and Q m is the matter density fraction. 
The cluster formation redshift zt is on average a decreasing function of the mass, 
meaning that smaller clusters are formed at larger redshifts, although fluctuations 
in the value of Zf from cluster to cluster at given mass are present due to the 
stochastic nature of the merger tree. 

The relative velocity of the two merging clusters, V r , can be easily calculated 
from energy conservation: 

Tvir,! + T v ir,2 * ^-(112 

where M r = M\Mij (M1+M2) is the reduced mass and R12 is the turnaround radius 
of the system, where the two subclusters are supposed to have zero relative velocity. 
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Fig. 2. LEFT PANEL: Merger history of a cluster with present mass 10 15 solar masses. The 
mass (y-axis) suffers major jumps in big merger events. Time is on the x-axis. RIGHT PANEL: 
Distribution of the Mach numbers of merger related shocks for 500 realizations of the merger tree 
of a cluster with present mass 10 15 Mq. 



In fact the final value of the relative velocity at the merger is quite insensitive to 
the exact initial condition of the two subclusters, being quite close to the free fall 
velocity. 

The sound speed in the cluster i is given by ^21 

2 t , GMi 

^'vir,i 

where the virial theorem has been used to relate the gas temperature to the mass 
and virial radius of the cluster. The adiabatic index of the gas is j g = 5/3. The 
Mach number of each cluster while moving in the volume of the other cluster is then 
Mi = V r /c s _i. Here it has been implicitely assumed that the Mach number remains 
constant during the whole merger. 

The procedure illustrated above can be applied to a generic couple of merging 
clusters, and in particular it can be applied to a generic merger event in the history 
of a cluster with fixed mass at the present time. 

In Fig. 2 (right panel) the differential distribution of shock Mach numbers is 
shown. Results have been extracte d fr om 500 realizations of the merger history of a 
cluster with present mass 10 15 .MqE2]. The most striking feature of the histogram is 
a very pronounced peak at Mach number ~ 1.4. This means that merger shocks are 
generally very weak, most of them are mildly supersonic, and only ~ 5% of merger 
shocks populate the tail of the distribution, with Mach number above ~ 3 ^1. 

The Mach number may be higher than that provided by the procedure illustrated 
above if the merger takes place in an overdense region such as a supercluster, or if 
the whole temporal evolution of the merger is considered instead of a mean value 
during the merger event. However, it has been demonstrated that none of these 
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effects alters significantly the results illustrated above: the former effect is negligible 
for very massive clusters ^§1 anc j the latter leads only to a slight increase of the 
merger velocity, less than a factor of two, but limited to a small fraction of the 
duration of the merger event ^31. Some additional uncertainties might arise from 
the assumption that the merging clusters are virialized to start with and from the 
assumptions made on the spatial distribution of dark matter in the parent halos. 

Due to the intrinsic complexity of the process of structure formation, numerical 
simulations are the best way to study properties of cosmological shock waves such 
as the shock velocity, strength and morphology, although much care is needed in 
the identification of shocks and in evaluating their efficiency as particle accelera- 
tors. Results from earlier simulations ^51 were in clear disagreement with the above 
mentioned semi-analytical findings. The distribution of Mach numbers reported in 
^Sl for shocks located in the inner 0.5h~ 1 Mpc region of clusters exhibited a pro- 
nounced peak at A4 ~ 4 -r 5, weaker shocks being virtually absent. More recently, 
higher resolution simulations led to a revision of these earlier results, and signif- 
icantly mitigated the discrepancy between the two approaches. According to 
there is a peak in the distribution of shock surfaces in the whole simulated box at 
A4 ~ 1.5, in very nice agreement with Fig. 2. Moreover, the same authors showed 
that most of the shock energy is dissipated at weak shocks (A4 ~ 2-^3). These 
results have been confirmed in^^, where energy is found to be mostly dissipated at 
shocks with Mach number ~ 2. 

X-ray observations have allowed the measurement of the Mach number of 
merger-related shocks in the inner regions of a few clusters. To date, the strongest 
shock ever observed in the ICM has a Mach number of the order of ~ 3 SSI, while 
claims of detection of weaker shocks have been reported l &7 | 88 | 47 [ 

Besides mergers, the growth of structures proceeds through the so-called sec- 
ondary infall onto already formed objects 1^21, or accretion. Clusters continuously 
attract unbound matter in their neighborhood and their mass increases due to ac- 
cretion of this new material. Accretion shocks form around clusters and propagate 
outwards, carrying the information on the virialization of the inner regions. The 
existence of non virialized gas outside clusters is suggested by cosmological simula- 
tions SSI. This gas is structured in low density, relatively cold filaments and sheets 
(T ~ 10 4 -j- 10 6 -?£T) whose X-ray emission is expected to be tenuous and hardly 
detectable by present day instruments. The low temperature of the gas implies that 
accre tion s hocks are always rather strong, with Mach numbers that may exceed 
~ 10 1 40 | 80| rpkg f rac tj on f the total cluster mass accreted through these shocks is 
less than ~ 10% ^2. This implies that, even if accretion shocks are very strong and 
thus able to dissipate efficiently the kinetic energy of the accret ion flow, they do 
not contribute significantly to the total energy budget of clusters [ 91 | 40 | 86 | 
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5. Diffusive particle acceleration at merger and accretion shocks 

A vast literatu re exists on th e problem of diffusive particle acceleration at colli- 
sionless shocks | 92 | 93 | 94 | 95 | 96 [ j n ^ g ^ es ^ p ar fo c i e version, the theory is simple and 
predicts that the spectrum of accelerated particles is a power law N(E) oc , 
with a slope which is determined uniquely by the Mach number M. of the supersonic 
fluid, 7 — 2^23^. For shocks with Mach number M. = 2 — 4 the spectrum of the 
accelerated particles varies between E~ z zz and E~ 2 - 27 . For M. = 1.4, where the 
peak of merger related shocks was found, the spectrum is as steep as E~ 62 : such 
weak shocks can only play a role as heaters of the intergalactic medium, while being 
irrelevant for particle acceleration. As we discuss below, the efficiency of injection at 
such weak shocks is also expected to be negligibly small. The convolution between 
the efficiency as particle accelerators and the steepeness of the spectrum selects the 
merger related shocks that are most relevant for particle acceleration in the ICM 

mi 

The process of acceleration of hadrons and electrons during cluster mergers has 
different peculiar aspects. The time scale for energy losses of electrons is shorter 
than the mean time between two mergers. Thus, it is a good approximation to 
assume that only the last merger event contributes to accelerate electrons and that 
these radiate for approximately their loss time. The spectrum of the accelerated 
electrons at each shock is a power law, with a slope which is expected to be hard 
enough to produce appreciable non-thermal emission only i n th e few (?s 5%) mergers 
characterized by a relatively high Mach number (A4 > 3) I^U 

The confinement of hadrons in the cluster volume for cosmological time scales 
implies that any new merger related shock front, besides accelerating a new pop- 
ulation of particles from the thermal pool, also re-energizes previously accelerated 
nuclei I^Sl. As a consequence, the spectrum of the hadrons present at any time in a 
cluster is the result of both acceleration and reacceleration of all previously confined 
hadrons. The combination of these effects makes the spectra of confined hadrons 
depart from a power law 1221, being steeper at low energies and gradually flatter 
at higher energies. This result was obtained by making the simple but most likely 
unrealistic assumption that the acceleration efficiency is independent of the Mach 
number (see below for a discussion of the problem based on the non-linear theory 
of shock acceleration). It is worth stressing that while the efficiency of a shock as a 
particle accelerator depends on the unknown details of injection, the efficiency for 
reacceleration is known I^ZI. 

The situation changes dramatically if one considers strong accretion shocks. 
According with the standard (test particle) theory of particle acceleration, these 
shocks should result in a universal spectrum cx E~ 2 since their Mach number is 
expected to be very large. On the other hand, in these conditions test-particle theory 
is known to fail, to an extent which depends upon several aspects of the problem 
(maximum energy, Mach number, presence or absence of turbulent heating, efficient 
generation of waves through streaming instability). The violations of test particle 
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Fig. 3. Proton spectra at the shock location for different values of the shock Mach number and 
of the maximum momentum of the accelerated particles. 

theory are grouped together in what is named shock modification, which is mainly 
determined by the dynamical reaction of the accelerated particles on the shock 
structure. 

The development of a non-linear theory of diffusive particle acceleration at shock 
fronts, which allows to accurately describe the shock modification and at the same 
time calculate the spectrum of the accelerated particles, is one of the major develop- 
men ts in the invest igation of the process of particle acceleration in the last decade 
(see | 92 | 93 | 94 | 95 | 9 6] f or re cent reviews). The main points of the non-linear theory 
of diffusive shock acceleration can be summarized as follows: 1) strong shocks are 
efficient accelerators, to the point that the dynamical reaction of the accelerated par- 
ticles cannot be neglected; 2) this dynamical reaction is responsible for a gradual 
slowing down of the fluid upstream of the shock front, the so-called shock precursor, 
which reflects the spatial distribution of the pressure of the accelerated particles; 
3) the precursor results in a non-power law, concave spectrum, flatter than the 
test-particle prediction (N(p) oc p~ 4 ) at high energies and steeper than that at low 
energies; 4) as a consequence of efficient particle acceleration, less energy is available 
for heating of the downstream gas, which therefore is predicted to be colder than 
expected on the basis of the standard Rankine-Hugoniot relations. 

The shock modification induced by the dynamical reaction of the accelerated 
particles leads to a determination of the efficiency for particle acceleration as a 
function of the Mach number of the shock, which is crucial for assessing the role of 
shock fronts for particle acceleration during structure formation. 

Some of the i mplications of the theory for large scale structures have been dis- 
cussed mJ 98 | 100 ] anc [ are now recognized as crucial to understand the role of cosmic 
rays in clusters of galaxies and their neighborhood. In order to illustrate the main 
results, we ad opt here the approach described in where injection was treated 
following Ref. 

In Fig. 3 we plot the spectra of protons at the shock location for Mach number 
M = 100 (upper panel) and M — 10 (lower panel) and for E max /(mc 2 ) = 10 7 and 
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Fig. 4. a) Temperature ratio between downstream and upstream infinity as a function of the Mach 
number (see text), b) Efficiencies for advected particles, escaping particles and total efficiency as 
functions of the Mach number. 

5 x 10 10 . These values of the Mach number are clearly suitable to describe accretion 
shocks rather than merger shocks. The maximum energy of particles accelerated 
at accretion shocks is believed to be within the range of values considered here 
| 101 | 102| 103| rjijjg normalization of the curves is such that these spectra reproduce 
those expected for protons accelerated at the accretion shock of a Coma-like cluster 
and therein confined ^1 The concave shape of the spectra, typical of non-linear 
particle acceleration at shocks, is clear in the figure. 

The Maxwellian distribution in Fig. 3 represents the spectrum of the particles in 
the downstream fluid, at the temperature which is calculated by solving the suitable 
conservation equations at the shock. As a comparison we also plot the prediction 
of test-particle theory assuming that a fraction of 10% of the shock kinetic energy 
flux is converted into cosmic rays (horizontal lines). 

As stressed above, one of the consequences of the non-linear backreaction of 
cosmic rays on the shocks induced by large scale structures is that the gas is heated 
less than it would be at shocks where cosmic rays are absent. In Fig. 4 (left panel), 
we plot the temperature ratio between downstream infinity and upstream infinity 
T2/T1 for an ordinary shock (dotted line) and for a cosmic ray modified shock, for 
E max = 10 2 , 10 6 and 5 x 10 10 GeV (from top to bottom) as a function of the Mach 
number of the shock. We can see that the temperature jump is basically unaffected 
by the presence of accelerated particles for low Mach numbers, but the heating is 
severely suppressed for strongly modified (high Mach number) shocks. 

In Fig. 4 (right panel) we also plot the fraction of flux in the form of relativistic 
particles (in units of (l/2)pu 3 ) which is effectively advected downstream, the frac- 
tion of flux which escapes the system at E max (from upstream) and the sum of the 
two, which saturates to a number very close to unity for large Mach numbers. Solid, 
dashed and dotted lines refer to E max = 5x 10 10 , 10 6 and 10 2 GeV respectively. 

It is worth stressing that the very poorly known phenomenon of turbulent heat- 
ing of the upstream plasma may somewhat reduce the shock modification with 
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respect to the results illustated above, though not affecting the basic conclusions. 

The non-linear effects of diffusive particle acceleration are know n to be crucial 
in the context of particle acceleration at supernova remnant shocks | 92 | 93 | 94 | 95 | 96 | 
In the case of clusters of galaxies these effects have never been explored in detail, 
all existing calculations being based on ver y simple assumptions on the spectrum 
of the accelerated part icles (see however ^Hl) , 

As discussed in 1^21 since non-linear effects in shock acceleration affect both the 
spectral shape and the overall normalization of the accelerated particles, modifi- 
cations in the predicted gamma ray fluxes from clusters are expected accordingly. 
This issue deserves further investigation, and it will likely constitute one of the most 
interesting future developments in the field. 



6. Gamma Rays from clusters 

The initial motivation for the interest in the gamma ray emission from clusters 

arose as a consequence of the possibility of cosmic ray confinement in the ICM, 

and consequently the possibility that radio halos could be the result of synchrotron 

emission from secondary electrons. A natural byproduct of the confinement is the 

copius emission of gamma radiation due to the production and decay of neutral 

pions. Calculations of this em ission from single clusters of galaxies have been carried 
out in 34I32I104I107I106TT051 

Gamma radiation in the ICM can also be generated as a result of ICS of high 
energy electrons off the universal photon background. In particular, electrons can 
be accelerated at cosmological shocks up to energies of tens of TeV. The resulting 
ICS emission extends to multi-TeV energies ^21. 

Other mechanisms of production of gamma rays in the ICM include ICS from 
secondary electrons, non-thermal bremsstrahlung and ICS from pairs generated in 
Bethe-Heitler processes between Ultra High Energy protons (if present) and photons 
in the cosmic microwave background. Th e fir st two channels may be relevant only 
for gamma ray energies below rj lOOMeV^ and will not be discussed further here. 
The latter mechanism will be briefly discussed in Sec. 16.31 

6.1. Gamma rays from tt° decay 

The present observational situation concerning radio halo emission from clusters is 
not favorable for secondary electron models, at least in those few cases in which the 
data are good enough to allow for such a type of analysis, as for the Coma cluster. In 
such cases the radio emission can be used to impose limits on the amount of hadronic 
cosmic rays trapped in the I CM and therefore on the flux of gamma rays due to tt° 
decay that may be expected Q23EI These limits are usually rather model dependent, 
in that the amount of cosmic ray energy that is allowed depends on magnetic field 
strength, slope of the proton spectrum and spatial distribution of cosmic rays in 
the intracluster medium. With these limitations, current limits suggest that not 
more than ~ 10 -r- 30% of the cluster thermal energy may be in form of relativistic 
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Fig. 5. a) Number of accreting (solid line) and merging (dashed line) clusters with gamma ray flux 
greater than F. The vertical lines represent the GLAST, AGILE and EGRET sensitivity for point 
sources, b) Gamma ray emission in the 100 GeV - 10 TeV region. The thick solid lines represent 
the sensitivities of a IACT for point sources (lower curve) and extended sources (upper curve). The 
predicted gamma ray fluxes from a Coma-like cluster at a distance of 100 Mpc with and without 
absorption of the infrared background are plotted as dashed and solid lines respectively. 

particles. If the amount of relativistic particles is not much smaller than this limit, 
GLAST s hould easily detect the gamma rays from pion decay from several rich 
clusters l 32 l 108 l in the energy range 100 MeV - 10 GeV. Other limits on the gamma 
ray emission of some nearby clu sters of galaxies were calculated by using the existing 
EGRET upper limits 1 104 | 1Q 7] 

Clusters of galaxies are however interesting targets even for Imaging Atmo- 
spheric Cherenkov Telescopes such as HESS, MAGIC and VERITAS. We briefly 
consider here the perspective for detection of a Coma-like cluster at TeV energies. 
If cosmic rays have a power law spectrum in the form N(E) = NqE~ 2 (for the pur- 
pose of estimating the flux we are neglecting here the non-linear effects discussed in 
the previous section), and are spatially distributed following the thermal gas, the 
expected gamma ray flux is E 2 F(E) » 10~ 12 (r]/0.1)erg/cm 2 /s, wher e 77 is the frac- 
tion of the cluster thermal energy in the form of relativistic protons [ 106 | 104 | The 
spectrum of the gamma radiation is also a power law at large enough energy, with 
a high energy cutoff which may be due to either lack of confinement of the parent 
protons or to absorption of gamma rays in the infrared background (in fact in the 
former case a steepening is expected rather than a cutoff). At the low redshift of 
the Coma cluster (z = 0.023) the o pacity due to absorption becomes of order unity 
at energies above several TeV 

The expected brightness profile should be centrally peaked, tracing the spatial 
distribution of the thermal gas, which represents the target for pp collisions. Thus, 
the emission is expected to be concentrated within a core of apparent size of ~ 
0.5°-f-l° for a Coma-like cluster. Taking as a reference the HESS sensitivity at 1 TeV 
for sources uniformly extended on i} ext ~ 1°, which is roughly ~ 3 x 10~ 12 erg/cm 2 /s 
for 10 hours exposure, we can estimate the on source exposure time needed for a 
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detection as a function of the parameter rj: 




Thus, a few tens of hours of exposure time seem adequate to detect a Coma-like 
cluster in the most optimistic case in which rj sa 0.1, while deeper observations 
(> 100 hours) are needed to probe more conservative flux levels (e.g. lower values 
of rj and/or a steeper proton spectrum). 

The detection of TeV photons with a brightness profile correlated with the gas 
density profile would hint to the presence of relativistic hadrons in the ICM and 
would allow to finally measure the cluster non thermal energy content. 



6.2. Gamma rays from Inverse Compton Scattering (ICS) 

The detectability of the ICS gamma ray signal from electrons accelerated at either 
merger or accretion related shocks crucially depends on the strength of the shocks 
formed in cither one of these cases. Only strong shocks can accelerate particles with 
spectra hard enough to result in an appreciable gamma ray emission. 

The LogN-LogS diagram of clusters of galaxies as gamma ray sources was calcu- 
lated misusing the PS based formalism described in Sec.|H The balance between 
injection and energy losses drives electrons towards their time independent equilib- 
rium distribution, with a spectrum one power steeper than the injection spectrum. 
The results of ^1 are shown in Fig. 5 (left panel), where the typical sensitivities 
of GLAST, AGILE and EGRET are also indicated. The calculations refer to the 
case of a constant efficiency of electron acceleration of 5%, independent of the Mach 
number of the shocks. It is worth stressing that the electrons accelerated according 
with this recipe provide a negligible contribution to the diffuse radio emission and 
to the hard X-ray emission. Within these assumptions, about ~ 50 clusters should 
be detected in gamma rays above 100 MeV by GLAST, equally distributed between 
merging and accreting clusters. AGILE on the other hand might be able to detect 
~ 10 clusters. Moreover, according to these predictions, no cluster should have been 
detected by EGRET, in agreement with what observed The self-consistent de- 
termination of the shock Mach numb ers and o f the slope of electron spectra is the 
crucial ingredient of these calculations ^ ^ I 83 J84] p rev ious approaches, relying on the 
incorrect assumption that all the structure formation shocks are strong (and thus 
the electron spectrum is always the canonical E~ 2 power law), led to exceedingly 
optimistic predictions on the number of detectable sources I 1H | H0 [ 

The efficiency of acceleration of electrons at shocks is an unknown quantity. The 
value of 5% which is generally adopte d is a reasonable (maybe even optimistic) 

1 1 1 01 

guess, but nothing more. Recent claims 1 1 that this efficiency could be measured 
in the case of supernova shocks also appears rather optimistic, given the numerous 
assumptions that such as estimate is based upon (for instance the Mach number 
of the shock is very different in the two cases, the shock modification could play a 
crucial role, the strength of the magnetic field could be dramatically different and 
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for remnants showing TeV emission it is hard to distinguish between leptonic and 
hadronic scenarios). 

In Fig. 5 (right panel) we show the ICS gamma ray spectrum in the TeV domain 
as expected from a Coma-like cluster of galaxies at 100 Mpc distance. The effect of 
the gamma ray absorption in the infrared background is illustrated by the difference 
between the solid lines (with absorption) and dashed lines (without absorption). 
From top to bottom, the lines refer to three different cases: 1) a merger between two 
clusters with masses 10 15 Mq and 10 13 M©; 2) an accreting cluster with mass 10 15 Afo 
with a magnetic field at the shock in the upstream region O.l^tG; 3) an accreting 
cluster with mass 10 15 Mq with a magnetic field at the shock in the upstream region 
0.01/iG. The thick solid lines r epres ent the sensitivities for a generic Cherenkov 
telescope array as calculated in ^3 f or a point like and a ~ 1° extended source, 
which is roughly the apparent size of two merging clusters at 100 Mpc distance. 
Accreting clusters are even more extended (~ 2 — 3°) and thus harder to detect. 
More specifically, a telescope like HESS could probe the flux levels considered here 
only after > 50 — 100 hours of observation. 



6.3. Gamma rays from interactions of Ultra High Energy protons 
in the ICM 

Another channel for the production of gamma rays in the ICM is related to the 
possi bility th at Ultra High Energy protons (E > 10 18 eV) may be present in the 
ICM H1 4 | 115 | Before escaping the cluster, a proton with energy E p loses a fraction 
of its energy by interacting with cosmic microwave background photons, producing 
electron-positron pairs via the Bethe-Heitler process. Pairs are produced with en- 
ergy ~ (2m e /m p )E p . These pairs rapidly cool via synchrotron emission and ICS, 
producing radiation in the hard-X and TeV gamma ray bands respectively. The 
spectrum of the radiation is a power law with a hard spectrum (photon index 
r ~ 1.5), so that no significant emission is expect ed in the GeV domain. Possible 
sources of very high energy protons include AGNs or cluster accretion shocks 
ESI The corresponding gamma ray flux is however extremely uncertain and in the 
absence of efficient confinement up to such high energies the flux may be too low 
to be ever detectable. 



7. Diffuse gamma ray background from structure formation 

In the last decade clusters of galaxies have often been proposed as the sources of the 
extragalactic diffuse gamma ray background, but more accurate calculations have 
always lowered their contribution to a negligible fractio n. This has been a hot to pic 
in the field, with numerous papers addressing the issue LLL 7 I 18 | 20 | 91 | 110 | 112 | 118 | i n 
this section we briefly summarize the relevant aspects of the problem. 

Early interest in clusters of galaxies as potential sources of the diffuse gamma 
ray background w as st im ulat ed by the discovery of cosmic ray confinement in the in- 
tracluster volume HHHH1 w hich increases the effective column density for the process 
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of pion production. In many potential sources of cosmic rays were considered, 
and the contribution of clusters of galaxies to the diffuse background observed by 
EGRET was foun d t o be negligible. This result was later confirmed in a more de- 
tailed calculation ^321, where however the spectrum was chosen ad hoc as a rather 
fiat power law (E~ a with a = 2.2 — 2.4), resembling the spectrum of particles ac- 
celerated at astrophysical shocks and having in mind the secondary electron mo del 
for the origin of radio halos. The diffuse flux of gamma rays as calculated in 
is plotted in Fig. [6] (Left panel), where the shaded region provides an estimate of 
the uncertainties arising from source luminosities and cosmological models. In the 
light of the most recent understanding of particle acceleration in clusters of galaxies, 
these fluxes should be considered as upper limits to the actual hadronic contribution 
of clusters of galaxies to the diffuse extragalactic gamma ray background. 




Fig. 6. Diffuse gamma ray background due to clusters of galaxies. Left Panel: contribution of 
hadronic interactions 32 with its uncertainty (shaded region). Right Panel: contribution from ICS 
of relativistic electrons. In both plots we also show the flux of diffuse extragalactic gamma ray 
background measured by EGRET 119 . 



As discussed in the previous sections, electrons are also expected to be accel- 
erated in the ICM. Here we concentrate our attention on electrons accelerated at 
either merger or accretion shocks. Electrons accelerated at each merger or accretion 
event lose energy through ICS on the photons of the cosmic microwave background, 
which get upscattered to gamma ray energies. 

The diffuse flux of gamma radiation due to ICS of electrons accelerated in merg- 
ers (dashed line) and during accretion (dotted line) is plotted in Fig. [6] (right panel) , 
where an acceleration efficiency (for electrons) of 5% h as be en assumed "^l The ob- 
served diffuse flux is represented by the shaded region ESI ( a recent re-analysis of 
the EGRET data led to a slightly lower extragalactic diffuse gamma ray background 
xhe calculation of the contribution of merger shocks requires to fix a minimum 
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mass of the merging clusters. The dashed line refers to a minimum mass ~ IO^Mq, 
while a more realistic result is plotted as a dash-dotted line, where the minimum 
mass is fixed at 1O 13 M0 (roughly the mass of a galaxy group). In the same figure we 
plot for comparison the predictions of^^ (solid line), where the same acceleration 
efficiency was adopted. Several arguments have been presented in to explain the 
difference between their predictions, illustrated in Fig. [6] and those of I 20|llt 3 ^ e 
basic points being, as already discussed above, that it is crucial to properly evaluate 
the Mach numbers of the merger shocks and that it is somewhat important to keep 
track of the redshift evolution of the merger rate. 

The contribution of merger shocks to the diffuse extragalactic gamma ray back- 
ground is therefore likely to be in between the dashed and the dash-dotted line. 

The flux of gamma radiation from both mergers and accretion is a factor ~ 
1 — 1 00 s maller than the observed background and smaller than the flux predicted 
i n | 20 | 110 ] k v th e same factor. An acceleration efficiency of the order of 50% should 
be adopted in order to reproduce observations. This would be unreasonable for 
electrons as accelerated particles, and would violate our initial assumption of shock 
acceleration in the linear regime (no backreaction of the accelerated particles on the 
shock). 

In a reevaluation of the diffuse gamma ray background from large scale 
structure shocks was carried out by using N-b ody simulations. The results seem to 
be consistent with those plotted in Fig.[6j InGSSthe authors emphasize the difficulties 
in the identification of shocks with Mach number below 10 (and the impossibility 
to detect shocks with Mach numbers below ~ 3 — 4) . 



8. Conclusions 

Clusters of galaxies are sites where non-thermal activity has been observed, mainly 
in the form of radio and hard X-ray radiation, resulting from synchrotron and in- 
verse Compton emission of relativistic electrons respectively. The sources of such 
electrons, the acceleration mechanisms involved in their production and the coexis- 
tence of a hadronic component in the ICM are subjects of intense debate, which we 
summarized here. At the present time observations seem to require a contribution 
from particle reacceleration through either resonant or non resonant particle-wave 
interactions. Although the general principle that this acceleration process is based 
upon is well understood, several details deserve more attention, especially in the 
direction of achieving a fully self-consistent picture of particle acceleration in the 
presence of plasma instabilities in the intracluster medium. The development of this 
line of investigation appears to be even more justified by the uprising of modern ob- 
servational facilities such as LOFAR and LWA, operating at low radio frequencies, 

II Ol I 

where giant radio halos are expected to be more frequent 1 I 

Both primary and secondary models of radio halos run into problems in explain- 
ing self-consistently morphological, spectral and statistical properties of extended 
radio halos. However, it should be clear that most of these constraints derive from 



February 5, 2008 15:34 



22 P. Blasi, S. Gabici & G. Brunetti 

direct analyses of the data on the Coma cluster and a few other clusters that have 
been studied in a systematic way, and might reveal themselves as incomplete when 
larger samples will become available. 

The detection of gamma rays from nearby clusters of galaxies would provide us 
with a powerful tool to discriminate among the various models proposed to explain 
observations. An overview of the existing predictions and models appears quite 
timely, since several gamma ray telescopes are starting to operate (HESS, MAGIC, 
VERITAS) and some others are coming soon (GLAST, AGILE). 

The basic mechanisms for the production of gamma radiation in the ICM are 
1) the decay of neutral pions produced in inelastic pp collisions and 2) ICS of high 
energy electrons on the photons of the cosmic microwave background. 

Both electrons and protons (or more in general nuclei) are expected to be ac- 
celerated at the shock fronts that form during the process of structure formation. 
These shocks are typically classified as merger shocks and accretion shocks, depend- 
ing on whether they are generated in the inner regions of the parent clusters (which 
makes them weaker due to the high temperature of the ICM) or in the outskirts 
and filaments, where fast, cold, non virialized plasma falls into the gravitational 
potential well of the parent cluster. A small fraction of merger shocks may also 
have relatively large Mach number. 

The contribution of the hadronic channel to gamma ray emission is rather model 
dependent because there may be several classes of sources that may contribute at 
different levels and with different spectra of accelerated particles. If acceleration 
takes place at the shocks that develop during structure formation, then an estimate 
of the energy in the form of non-thermal particles may be obtained from the merger 
history of a cluster as described by the hierarchical model of structure formation. 
At each merger event that a cluster suffers a new population of hadrons is i njected 
but the pre-existing cosmic ray population, diffusively confined in the ICM | 18 | 19 | 
is also reacceleratedlSH. This combination of effects' was shown to result in non- 
power law spectra. Due to the weakness of shocks in the inner parts of galaxy 
clusters, the spectra of accelerated particles are typically steep, and the proton 
spectrum accumulate d in th e ICM as a result of mergers does not lead to intense 
gamma ray emission I 39 | 84 [ Acceleration at accretion shocks in the outskirts of 
clusters is more promising (the Mach number is ~ 10 — 100) and the particles 
could be advected into the cluster and trapped there for cosmological times. The 
gamma ray emission due to ir° decay in this scenario might be detected by GLAST, 
depending on the acceleration efficiency at the accretion shock. If additional hadrons 
are injected by sources in the ICM (galaxies, AGNs) the gamma ray flux due to 
hadronic interactions may be somewhat higher tL 9 | 32 | 34 | p^^^gj. strong limits on the 
energy content allowed to be in the form of non-thermal hadrons can be imposed 
(mostly ado pting the Coma cluster as the prototype cluster) by using the EGRET 
upper limits I 32 | 104 | 37 | 107 |^ U pp er limits on the flux of TeV photo ns and the 

■ i noi 

existence of the cutoff in the radio spectrum at GHz frequency 1 I Typically these 
methods lead to conclude that not more than ~ 10 — 30% of the energy of the 
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cluster can be in the form of a non-thermal population of hadrons. These limits are 
however rather model dependent and should be used with care. If the amount of 
hadrons is not much smaller than this value, and particle spectra are hard enough 
(oc E~ 2 ), forthcoming observations in the TeV domain might finally detect the first 
gamma rays from clusters. 

The detection of gamma rays generated through ICS of high energy electrons 
accelerated at merger and accretion shocks appears somewhat more promising. As 
illustrated in Fig. [6j ~ 10 — 50 clusters should be detectable with GLAST in the 
energy range > 100 MeV. T his pr ediction is consistent with the fact that EGRET 
has not detected any cluster I 83 | 17 [ rp ne detection of gamma rays in the TeV energy 
range appears more problematic, due to the extended size of the emission region. 

The contribution of clusters of galaxies to the extragalactic diffuse gamma ray 
background is less than 1 — 10% of the flux measured in ^^ l As better observations 
will become available, it is possible that the classes of sources that contribute to 
the extragalactic gamma ray background may become better identified, either by 
subtraction of known sources, or by identification of new objects with a given gamma 
ray emission. 

From the theoretical point of view the most important development in the un- 
derstanding of the non-thermal aspects of structure formation will most likely come 
from the introduction of non linear effects in diffusive particle acceleration at shocks. 
These effects are showing all their importance for the case of supernova remnants 
as cosmic ray accelerators and it is likely that a similar role will arise in the future 
for large scale structure shocks. 

Acknowledgements 

The authors are very grateful to F. Aharonian, T. Jones, O. Reimer and E. Resconi 
for constructive comments on the manuscript. SG gratefully acknowledges support 
from the Alexander von Humboldt foundation. The research work of PB and GB 
was partially funded through grant PRIN2004. GB acknowledges partial support 
from grants PRIN2005 and PRIN-INAF2005. 

References 

1. C. L. Sarazin, X-ray emission from clusters of galaxies, Cambridge: Cambridge Uni- 
versity Press (1988). 

2. L. Feretti, Adv. Space Res. 36, 729 (2005). 

3. S. Giacintucci, T. Venturi, G. Brunetti, S. Bardelli, D. Dallacasa, S. Ettori, A. 
Finoguenov, A. P. Rao, and E. Zucca, A&A 440, 867 (2005) 

4. R. Fusco-Femiano, Astrophys. Space Set. 294, 37 (2004). 

5. R. Lieu, and J. P. D. Mittaz, in N. J. C. Spooner and V. Kudryavtsev eds., The 
Identification of Dark Matter, Singapore: World Scientific Publishing (2005). 

6. S. Bowyer, J. Korean Astron. Soc. 37, 295 (2004). 

7. G. Giovannini, and L. Feretti, in L. Feretti, I. M. Gioia and G. Giovannini eds., Merging 
Processes in Galaxy Clusters, Kluwer Academic Publishers, Dordrecht (2002). 



February 5, 2008 15:34 



24 P. Blast, S. Gabici & G. Brunetti 

8. G. Giovannini, M. Tordi, and L. Feretti, New Astron. 4, 141 (1999). 

9. P.-H. Kuo, C.-Y. Hwang, and W.-H. Ip, Astrophys. J. 604, 108 (2004). 

10. V. Petrosian, Astrophys. J. 557, 560 (2001). 

11. C. L. Carilli, and G. B. Taylor, Ann. Rev. Astron. Astrophys. 40, 319 (2002). 

12. F. Govoni, and L. Feretti, Int. J. Mod. Phys. D 13, 1549 (2004). 

13. T. E. Clarke, P. P. Kronberg, and H. Bohringer, Astrophys. J. 547, Llll (2001). 

14. R. Fusco-Femiano, M. Orlandini, G. Brunetti, L. Feretti, G. Giovannini, P. Grandi, and 
G. Setti, Astrophys. J. 602, L73 (2004); Y. Rephaeli, D. Gruber, Y. Arieli, Astrophys. 
J. 649, 673 (2006). 

15. W. I. Newman, A. L. Newman, and Y. Rephaeli, Astrophys. J. 575, 755 (2002). 

16. J. Bagchi, F. Durret, G. B. Lima-Neto, and S. Paul, Science 314, 791 (2006). 

17. O. Reimer, M. Pohl, P. Sreekumar, and J.R. Mattox, Astrophys. J. 588, 155 (2003). 

18. V. S. Berezinsky, P. Blasi, and V. S. Ptuskin, Astrophys. J. 487, 529 (1997). 

19. H. J. Volk, F. A. Aharonian, and D. Breitschwerdt, Space Set. Rev. 75, 279 (1996). 

20. A. Loeb, and E. Waxman, Nature 405 156 (2000). 

21. D. G. Wentzel, Ann. Rev. Astron. Astrophys. 12, 71 (1974). 

22. F. Govoni, Astron. Nach. 327, 539 (2006). 

23. C. Vogt, T. A. Ensslin, Astron. Astrophys. 434, 67 (2005). 

24. M. Murgia, F. Govoni, L. Feretti, G. Giovannini, D. Dallacasa, R. Fanti, G. B. Taylor, 
K. Dolag, Astron. Astrophys. 424, 429 (2004). 

25. W. J. Jaffe, Astrophys. J. 212, 1 (1977). 

26. G. Brunetti, JKAS 37, 493 (2004). 

27. B. Dennison, Astrophys. J. 239, L93 (1980). 

28. P. Blasi, and S. Colafrancesco, Astropart. Phys. 12, 169 (1999). 

29. R. Schlickeiser, A. Sievers, and H. Thiemann, Astron. Astrophys. 182, 21 (1987). 

30. P. C. Tribble, MNRAS 263, 31 (1993). 

31. G. Brunetti, G. Setti, L. Feretti, and G. Giovannini, MNRAS 320, 365 (2001). 

32. S. Colafrancesco, and P. Blasi, Astropart. Phys. 9, 227 (1998). 

33. H. J. Volk, and A. M. Atoyan, Astropart. Phys. 11, 73 (1999). 

34. T. A. Ensslin, P. L. Biermann, P. P. Kronberg, and X.-P. Wu, Astrophys. J. 477, 560 
(1997). 

35. R. Cen, and J. P. Ostriker, Astrophys. J. 514, 1 (1999). 

36. F. Miniati, D. Ryu, H. Kang, and T. W. Jones, Astrophys. J. 559, 59 (2001). 

37. P. Blasi, Astropart. Phys. 15, 223 (2001). 

38. C. L. Sarazin, in L. Feretti, I.M. Gioia and G. Giovannini eds., Merging Processes in 
Galaxy Clusters, Kluwer Academic Publishers, Dordrecht (2002). 

39. S. Gabici, and P. Blasi, Astrophys. J. 583, 695 (2003). 

40. D. Ryu, H. Kang, E. Hallman, and T. W. Jones, Astrophys. J. 593, 599 (2003). 

41. R. Cassano, and G. Brunetti, MNRAS 357, 1313 (2005). 

42. Y. Fujita, and C. L. Sarazin, Astrophys. J. 563, 660 (2001). 

43. M. Takizawa, and T. Naito, Astrophys. J. 535, 586 (2000). 

44. M. Markevitch, A. Vikhlinin, Astrophys, J. 563, 95 (2001). 

45. M. Markevitch, A. H. Gonzales, L. David, A. Vikhlinin, S. Murray, W. Forman,C. 
Jones, and W. Tucker, Astrophys. J. 567, L27 (2002). 

46. K. Rottiger, J. O. Burns, and J. M. Stone, Astrophys, J. 518, 603 (1999). 

47. M. Markevitch, F. Govoni, G. Brunetti, D. Jerius, Astrophys, J. 627, 733 (2005). 

48. T. E. Clarke, T. A. Ensslin, Astron. J. 131, 2900 (2006). 

49. E. Bertschinger, Astrophys. J. Suppl. 58, 39 (1985). 

50. F. Miniati, T. W. Jones, D. Ryu, and H. Kang, Astrophys. J. 562, 233 (2001). 

51. T. A. Ensslin, Copal-Krishna, Astron. Astrophys. 366, 26 (2001). 



February 5, 2008 15:34 



Gamma rays from clusters of galaxies 25 

52. M. Hoeft, M. Briiggen, G. Yepes, MNRAS 347, 389 (2004). 

53. K. Dolag, and T. A. Ensslin, Astron. Astrophys. 362, 151 (2000). 

54. D. Buote, Astrophys. J. 553, L15 (2001). 

55. P. Schuecker, H. Bohringer, T. H. Reiprich, and L. Feretti, Astron. Astrophys. 378, 
408 (2001). 

56. M. Thierbach, U. Klein, and R. Wielebinski, Astron. Astrophys. 397, 53 (2003). 

57. G. Giovannini, L. Feretti, T. Venturi, K.-T. Kim, P. P. Kronberg, Astrophys, J. 406, 
399 (1993). 

58. L. Feretti, E. Orrii, G. Brunetti, G. Giovannini, N. Kassim, G. Setti, Astron. Astro- 
phys. 423, 111 (2004). 

59. P. Blasi, MAS 37, 483 (2004). 

60. E. Fermi, Phys. Rev. 75, 1169 (1949). 

61. E. Fermi, Astrophys. J. 119, 1 (1954). 

62. G. Brunetti, and P. Blasi, MNRAS 363, 117 (2005). 

63. H. Ohno, M. Takizawa, and S. Shibata, Astrophys. J. 577, 658 (2002). 

64. Y. Fujita, M. Takizawa, and C. L. Sarazin, Astrophys. J. 584, 190 (2003). 

65. G. Brunetti, P. Blasi, R. Cassano, and S. Gabici, MNRAS 350, 1174 (2004). 

66. P. Goldreich, and S. Sridhar, Astrophys. J. 438, 763 (1995). 

67. J. Cho, and A. Lazarian, MNRAS 345, 325 (2003). 

68. H. Yan, and A. Lazarian, Astrophys. J. 614, 757 (2004). 

69. S. Kato, PAS' J 20, 59 (1968). 

70. A. Beresnyak, and A. Lazarian, Astrophys. J. 640, L175 (2006). 

71. J. Cho, and A. Lazarian, Astrophys. J. 638, 811 (2006). 

72. G. Brunetti, and A. Lazarian, submitted to MNRAS. 

73. A. Lazarian, Astrophys. J. 645, L25 (2006). 

74. G. Brunetti, Astron. Nachr. 327, 615 (2006). 

75. D. A. Buote, in L. Feretti, I. M. Gioia and G. Giovannini eds., Merging Processes in 
Galaxy Clusters, Kluwer Academic Publishers, Dordrecht (2002). 

76. W. Forman, C. Jones, M. Markevitch, A. Vikhlinin, in Merging Processes in Galaxy 
Clusters, eds. L. Feretti, I. M. Gioia and G. Giovannini (Kluwer Academic Publishers, 
Dordrecht, 2002) 

77. W. H. Press, and P. Schechter, Astrophys. J. 187, 425 (1974). 

78. J. R. Bond, S. Cole, G. Efstathiou, and N. Kaiser, Astrophys. J. 379, 440 (1991). 

79. C. Lacey, and S. Cole, MNRAS 262, 627 (1993). 

80. A. Cavaliere, N. Menci, and P. Tozzi, MNRAS 308, 599 (1999). 

81. S. J. Randall, C. L. Sarazin, and P. M. Ricker, Astrophys. J. 577, 579 (2002). 

82. M. Takizawa, Astrophys. J. 520, 514 (1999). 

83. S. Gabici, and P. Blasi, Astropart. Phys. 20, 579 (2004). 

84. R. C. Berrington, and C. D. Dermer, Astrophys. J 594, 709 (2003). 

85. F. Miniati, D. Ryu, H. Kang, T. W. Jones, R. Cen, and J. P. Ostriker, Astrophys. J. 
542, 608 (2000). 

86. C. Pfrommer, V. Springel, T. A. Ensslin, and M. Jubelgas, MNRAS 367, 113 (2006). 

87. M. Markevitch, C. L. Sarazin, and A. Vikhlinin, Astrophys. J. 521, 526 (1999). 

88. R. A. Krivonos, A. A. Vikhlinin, M. L. Markevitch, and M. N. Pavlinsky, Astron. Lett. 
29, 425 (2003). 

89. V. Pavlidou, and B. D. Fields, Astrophys. J. 642, 734 (2006). 

90. U. Keshet, E. Waxman, and A. Loeb, Astrophys. J. 617, 281 (2004). 

91. S. Gabici, and P. Blasi, Astropart. Phys. 19, 679 (2003). 

92. L. O'C. Drury, Rep. Prog. Phys. 46, 973 (1983). 

93. R. Blandford, and D. Eichler, Phys. Rep. 154, 1 (1987). 



February 5, 2008 15:34 



26 P. Blast, S. Gabici & G. Brunetti 

94. F. C. Jones, and D. C. Ellison, Space Set. Rev. 58, 259 (1991). 

95. M. A. Malkov, and L. O'C. Drury, Rep. Prog. Phys. 64, 429 (2001). 

96. P. Blasi, Mod. Phys. Lett. A20, 3055 (2005). 

97. A. R. Bell, MNRAS 182, 443 (1978). 

98. S. Gabici, and P. Blasi, in F. A. Aharonian, H. J. Volk and D. Horns eds., High 
Energy Gamma-Ray Astronomy, AIP Conference Proceedings, Volume 745. New York: 
American Institute of Physics, 2005., p. 561-566 

99. P. Blasi, S. Gabici, and G. Vannoni, MNRAS 361, 907 (2005). 

100. H. Kang, and T. W. Jones, Astrophys. J. 620, 44 (2005). 

101. H. Kang, J. P. Rachen, and P. L. Biermann, MNRAS 286, 257 (1997). 

102. H. Kang, D. Ryu, and T. W. Jones, Astrophys. J. 456, 422 (1996). 

103. M. Ostrowski, and G. Siemieniec-Ozieblo, Astron. Astrophys. 386, 829 (2002). 

104. P. Blasi, Astrophys. J. 525, 603 (1999). 

105. F. Miniati, MNRAS 342, 1009 (2003). 

106. A. M. Atoyan, and H. J. Volk, Astrophys. J. 535, 45 (2000). 

107. C. Pfrommer, and T. A. Ensslin, Astron. Astrophys. 413, 17 (2004). 

108. A. Reimer, O. Reimer, R. Schlickeiser, and A. Iyudin, Astron. Astrophys. 424, 773 
(2004). 

109. J. R. Primack, J. S. Bullock, and R. S. Somerville, in F. A. Aharonian, H. J. Volk 
and D. Horns eds., High Energy Gamma-Ray Astronomy, AIP Conference Proceedings, 
Volume 745. New York: American Institute of Physics, 2005., p.23-33. 

110. T. Totani, and T. Kitayama, Astrophys. J. 545, 572 (2000). 

111. E. Waxman, and A. Loeb, Astrophys. J. 545, Lll (2000). 

112. U. Keshet, E. Waxman, A. Loeb, V. Springel, and L. Hernquist, Astrophys. J. 585, 
128 (2003). 

113. F.A. Aharonian, W. Hofmann, A.K. Konopelko, and H.J. Volk, Astropart. Phys. 6, 
369 (1997). 

114. F. A. Aharonian, MNRAS 332, 215 (2002). 

115. C. Rordorf, D. Grasso, and K. Dolag, Astropart. Phys. 22, 167 (2004). 

116. S. Inouc, F. A. Aharonian, and N. Sugyiama, Astrophys. J. 628, L91 (2005). 

117. A. Dar, and N. J. Shaviv, Phys. Rev. Lett. 75, 3052 (1995) 

118. F. Miniati, MNRAS 337, 199 (2002). 

119. P. Sreekumar et al., Astrophys. J. 494, 523 (1998). 

120. A. W. Strong, I. V. Moskalenko, and O. Reimer, Astrophys. J. 613, 956 (2004). 

121. R. Cassano, G. Brunetti, and G. Setti, MNRAS 369, 1577 (2006). 



